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On the dynamic efficiency of internal shocks in magnetized 
relativistic outflows 
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ABSTRACT 

We study the dynamic efficiency of conversion of kinetic-to-thermal/magnetic energy of in- 
ternal shocks in relativistic magnetized outflows. We model internal shocks as being caused 
by collisions of shells of plasma with the same energy flux and a non-zero relative velocity. 
The contact surface, where the interaction between the shells takes place, can break up either 
into two oppositely moving shocks (in the frame where the contact surface is at rest), or into 
a reverse shock and a forward rarefaction. We find that for moderately magnetized shocks 
(magnetization cr 0.1), the dynamic efficiency in a single two-shell interaction can be as 
large as 40%. Thus, the dynamic efficiency of moderately magnetized shocks is larger than in 
the corresponding unmagnetized two-shell interaction. If the slower shell propagates with a 
sufficiently large velocity, the efficiency is only weakly dependent on its Lorentz factor. Con- 
sequently, the dynamic efficiency of shell interactions in the magnetized flow of blazars and 
gamma-ray bursts is effectively the same. These results are quantitatively rather independent 
on the equation of state of the plasma. The radiative efficiency of the process is expected to 
be a fraction f r < 1 of the estimated dynamic one, the exact value of f r depending on the 
particularities of the emission processes which radiate away the thermal or magnetic energy 
of the shocked states. 

Key words: Hydrodynamics - (magnetohydrodynamics) MHD - Shock waves - gamma- 
rays: bursts - galaxies: BL Lacertae objects: general 



1 INTRODUCTION 



Internal shocks jRees & MeszarosI 1 19941) are invoked to explain 
the va riability of blazars (see e.g., Spada et al]|200ll : IMimica et al.l 



1 20051) and the light curves of the prompt phase of gamma-ray 
bursts (GRBs) jSari & Pirar]|l995l. Il997r . IDaigne & Mochkovitchl 
1 19981) . A possible problem in this model is the question whether 
this mechanism is efficient enough to explain the relation between 
the observed energi es both in the prompt GRB pha se and in the 
afterglow (see e.g.. IKobavashi et al.ll 19971 (KPS9 7), iBeloborodovl 
12000. IKobavashi & SarfeOOOllFan & Pirarjl 20061) . To asses the ef- 
ficiency of the internal shock model, most of the previous works 
focus on the comparison between the observed light curves and the 
model predictions emp loying a simple ine l astic collision of two- 
point masses (KPS97: ILazzati et all 1 19991; iNakar & Pirarj 120021 : 
iTanihata et a l. 2003; Zhan g & Meszarosll2004l) . Less attention has 
bee n paid to the hydrodynamic effects during the shell collision (but 
see IKobavashi & Sarill2000l : iKino et alj|200i iMimica et al1l2004 
l2005l : lBosniak et al^OCm 

The ejecta in GRBs and blazars may be rather magnetized, 
particularl y if they ar e originated as a Poynting-fiux-dominated 
flow (e.g., lusov|[T992h Forming shocks in highly magnetized me- 
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dia is challenging (Ree s&Gunrilll974l: l Kennel & Coronitil [l984T) . 
Therefore, to account for the observed phenomenology it is nec- 
essary to address how efficient the process of internal collisions in 
arbitrarily magnetized flows is. This question has been partly con- 
sidere d by a few recent works (e.g. jFan et al.|[2004l : IMimica et al.l 
120071) . 

The base to study the efficiency of internal collisions is the 
determination of the dynamic efficiency of a single binary collision, 
i.e., the efficiency of converting the kinetic energy of the colliding 
fluid into thermal and/or magnetic energy. Note that the radiative 
efficiency (i.e., the efficiency of converting the kinetic energy of the 
flow int o radiation) is expected t o be somewh a t smal ler. According 
to, e.g., IPanaitescu et al .1 J 1 9991) and iKumad J 1 999h . it can be as 
low as f r ~ 0.1. As we shall show in this paper, binary collisions 
in relativistic, magnetized flows can be an efficient enough way to 
dissipate a major fraction of the bulk kinetic energy of a relativistic 
ejecta. Therefore, it will depend on the efficiency of the particular 
radiation mechanism, that produces the observed emission (i.e., on 
the factor /,), that the model of internal shocks be efficient enough 
to explain the observations (particularly, the distribution of energies 
between the prompt GRB phase and the afterglow phase). 

We model internal shocks as shells of plasma with the same 
energy flux and a non-zero relative velocity. The contact surface, 
where the interaction between the shells takes place, can break up 
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either into two oppositely moving shocks (in the frame where the 
contact surface is at rest), or into a reverse shock and a forward rar- 
efaction. The determination of whether one or the other possibil- 
ity occurs is computed by estimating the invariant relative velocity 
between the fastest and the slowest shell, i.e., by solving the Rie- 
mann problem posed by the piecewise uniform states given by the 
physical quantities on the two interacting shells (Section[2]l. In Sec- 
tion[3]we define precisely the notion of dynamic efficiency, both for 
shocks and rarefactions. We perform a parametric study of the bi- 
nary shell collision dynamic efficiency in Section|4] The discussion 
and conclusions are listed in Section [5] Finally, we have extended 
our analysis on the dynamic efficiency of internal shocks in mag- 
netized, relativistic plasma to consider more realistic equations of 
state in the Appendix. 



2 RELATIVISTIC MAGNETOHYDRODYNAMIC 
RIEMANN PROBLEM 

We model the interaction between parts of the outflow with varying 
properties by considering Riemann problems, i.e. relativistic mag- 
netohydrodynamic initial-value problems with two constant states 
separated by a discontinuity in planar symmetry. We note, that we 
could use a more sophisticated approach consisting on performing 
numerical relativistic magnetohydrodynamic (RMHD) simulations 
of the interaction of parts of the outflow with different velocities. 
However, such an approach demands huge computational resources 
(eve n performing one dim ensional simulations using the same code 
as in iMimica et ai]|2009r) . and we are interested in sampling very 
finely a large parameter space with our models. Apart from this nu- 
merical reason, it is in order to point out that, by the internal shock 
phase, the lateral expansion of the flow is very small, since the flow 
is probably cold and ultrarelativistic. Thus, a description of the in- 
teractions assuming planar symmetry suffices to compute the dy- 
namic efficiency of such interactions (rather than a more complex 
spherically symmetric approach). 

In the following we use subscripts L and R to denote prop- 
erties of the (faster) left and (slower) right state, respectively. To 
avoid repeated writing of a factor 4n and the speed of light c, we 
normalize the rest-mass density p to p R , the energy density to p R c 2 
and the magnetic field strength to c ^jAnp R . 



2.1 Initial states of the Riemann problem 

For the initial thermal pressure of both states we assume that it 
is small fraction of the density, p L = xPl and p R = x- We as- 
sume magnetic fields perpendicular to the direction of the flow 
propagation. The remaining parameters determining the RMHD 
Riemann problem are: the density contrast p L , the Lorentz fac- 
tor of the right state T R , the relative Lorentz factor difference 
Ag := (T L - T R )/T R , and the magnetizations of left and right states, 
Ol := B 2 J(J%{\ + Ag) 2 p L ) and cr R := B 2 R /r 2 R , where B L and B R 
are the lab frame magnetic field strengths of left and right states, 
respectively. Furthermore, we define the total (thermal + magnetic) 
pressure 



■p(l + e) + 



ap 



(3) 



ap 

the total specific enthalpy 
h* := 1 + e + pip + a , 
and 



(1) 



(2) 



where e denotes the specific internal energy and is dependant on 
the equation of state used (see Section lTTI i. 

The general solution of a RMH D Riemann problem was found 
by iGiacomazzo & Rezzollal d2006t). and recently used in RMHD 
numerical codes by e.g. . Ivan der Hoist et al. I d2008l) . However, here 
we deal with a de generate RMHD con figuration, which solution 
was first found by iRomero et al.l d2005l) . The typical structure of 
the flow after the break up of the initial discontinuity consists of 
two initial states, and two intermediate states separated by a contact 
discontinuity (CD). The total pressure and velocity are the same on 
both sides of the CD. The quantity ajp is uniform everywhere, ex- 
cept across the CD, where it can have a jump. We denote the total 
pressure of intermediate states p" s , and rest-mass density left and 
right of the CD as ps_ L and ps tR . In the context of internal shocks, if 
the flow is ultrarelativistic in the direction of propagation, the ve- 
locity components perpendicular to the flow propagation should be 
negligibly small and, hence, they are set up to zero in our modeQ. 



2.2 Conditions for the existence of a two-shock solution 

One of the key steps in solving a Riemann problem is to determine 
under which conditions internal shocks can form. States ahead and 
behind the shock fr ont are related by the Lichnerowicz adiabat 
dRomero et alj2005b 

(hi) 2 - (K) 2 - I -i. - M (p b -p' a ) = 0. (4) 

\Pb Pa ) 

Following iRezzolla & Zanottil d200ll) . we study the relative veloc- 
ity between the states ahead (a) and behind (b) the shock front (all 
velocities are measured in the rest frame of the shock, and all ther- 
modynamic properties are measured in the fluid rest frame), 



Vab 



1 - V a Vj, 



(p* b - pPK - O 

(e* + p* b )(e* b + p*) 



(5) 



In our case states ahead of the shock are the initial (L, R) 
states, while states behind the shock are the intermediate states. 
Since v„i, is Lorentz-invariant, we can measure the velocity ahead 
of the left-propagating {reverse) shock (RS) in the frame in which 
the CD is at rest, 



Vi 



(P' S -pl)( e 's,L(p' S )- e l) 
(el+p* s )(e* SL (j>* s ) + pl) 



(6) 



Likewise, the velocity ahead of the right-going (forward) shock 
(FS) measured in the CD frame is 



(Ps -PrMApV-Q 
(e R + p* s )(e* (p* s ) + p R ) ' 



(7) 



where e* L and et R are the energy densities of the states to the left 
and to the right of the CD, respectively. The rest-mass densities p SJ1 
and p L R can be obtained from I0 and {2). 

Since both FS and RS only exist if pi > p* and p* s > pi, 
respectively, with decreasing p* s either the FS will disappear first 
(for pj = p* R > p* L , giving v r = 0) or the RS will disappear first (for 
Ps = Pl > P*R' gi ym g v i = 0)- Using equations ® and l[7} and the 

' If such velocities were significant, appre ciable changes in the R ie- 
mann structure may r esult as pointed out in|Alov & Rezzolkl i200d) or 
lAlov & Mimicl feOOgl) . 
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invariance of the relative velocity between the left and right states, 
v ir '■= (v; - v r )/(l - V;V r ), we can determine the minimum relative 
velocity for which a two-shock solution is possible 



(Vlr)lS 



(Pi - Pr^'sJpD ~ e P 
(e* SJt (PL*) + p' R )(e' R + p'j) 

(P' R - P"i)(el l(Pr) - e D 
(e s Ap R ) + Pl){e' L + p' R ) 



if pi = p* s > p' R 



if P'l<P'r = P*s 



(8) 



Generally, the quantity (v ; ,.) 2 s can t> e on ly determined numer- 
ically. If (v; r ) < (v lr )2s , a single shock and a rarefaction emerge 
from the initial discontinuity. It is even possible that instead of two 
shocks two rarefactions form (see, Rezzolla & Zanotti 200 lh . 



3 ENERGY DISSIPATION EFFICIENCY OF INTERNAL 
SHOCKS 

Internal shocks in relativistic outflows are invoked as moderately 
efficient means of conversion of the kinetic energy of the flow into 
radiation. In this section we present our model for inhomogeneous, 
ultrarelativistic outflows and provide an operative definition for the 
efficiency of conversion of the initial energy of the outflow into 
thermal and magnetic energy produced by internal shocks. We as- 
sume that a fraction of this thermal and magnetic energy will be 
radiated away. 

3.1 Outflow model 

To study internal shocks we idealize interactions of parts of the out- 
flow moving with different velocities as collisions of homogeneous 
shells. In our model the faster (left) shell catches up with the slower 
(right) one yielding, in some cases, a pair of shocks propagating 
in opposite directions (as seen from the CD frame). In order to 
cover a wide range of possible flow Lorentz factors and shell mag- 
netizations, we assume that initially, the flux of energy in the lab 
frame is uniform and the same in both shells! The energy flux for 
a shell with rest-mass density p, ratio of thermal pr essure to density 
X, ma gnetization a and Lorentz factor T is (e.g.. iLeismann et"aD 
I2005T) . 

f t := P [r 2 (i + e+^ + o-)-r] Vi-r- 2 . (9) 

Using the notation introduced in Section 12.11 and assuming the 
equality of F T in both shells we find that the density contrast p L 
between left and right shells is 



PL 



d + Asr^i + f+^ + ^-iV] ^i - r- 2 

[l + e+x + a-L- I?(l + Ag)- 1 ] ^/l-r^O+Ag)- 



(10) 



Considering cr L , cr R , T R and Ag as parameters, we can use dlOt 
to compute the rest of the variables needed to set-up the Riemann 
problem. We then compute the break up of the initial discontinuity 
between both shells using the exact Riemann solver developed by 
iRomero et aIU2005l) . 

In the following we use a polytropic equation of state with an 
adiabatic index y = 4/3: 



2 In|Mimica et al. <2009h we use a similar model to compare afterglow 
ejecta shells with different levels of magnetization, with a slight difference 
that in that study, instead of having the same flux of energy, all the ejecta 
shells were assumed to contain the same total energy. 



(r - i)p 

As we show in the Appendix, the Riemann solver of lRomero et al .1 
d2005h has been suitably modified to include a more realistic equa- 
tion of state. However, we find no qualitative differences between 
the results using a polytropic EoS (with either y = 4/3ory = 5/3) 
and y- variable EoS. Furthermore, the quantitative differences are 
very small in terms of dynamic efficiency. 

3.2 Efficiency of energy dissipation by a shock 

To model the dynamic e fficiency of energy dis sipation we follow 
the approach described in lMimica et al .1 1 2007t) . suitably modified 
to account for the fact that in the present work, there can occur sit- 
uations where either the FS or the RS do not exist (see Section l2~2l l. 
By using the exact solver we determine the existence of shocks and 
(in case one or two shocks exist) obtain the hydrodynamic state 
of the shocked fluid. We use subscripts S,L and S,R to denote 
shocked portions of left and right shells, respectively. In the fol- 
lowing we treat the efficiency of each shock separately. 



3.2.1 Reverse shock 

To compute the efficiency we need to compare the energy content 
of the initial (unshocked) faster shell with that of the shocked shell 
at the moment when RS has crossed the initial shell. Assuming an 
initial shell width Ax, we define total initial kinetic , thermal and 
magn etic energy (see also equations (A. 1) - (A.3) of lMimica et al .1 
120071) 



E K (T,p,Ax) 



E T (T,p,p, Ax) 



E M (T,p,o-,Ax) 



T(r - \)pAx 



[{pe + p)T 2 - p]Ax 



r 2 per Ax 



(12) 



When the RS crosses the whole initial shell, the length of the com- 
pressed shell (i.e., the fluid between the RS and the CD) is ( L Ax, 
where 

VCD ~ Vs,L . 

(l ■= < 1 

VL - V S ,L 

and x'cd and v$,l are velocities (in the lab frame) of the contact 
discontinuity and the RS, both obtained from the solver. Without 
loss of generality, we can normalize the initial shell width so that 
Ax = I. Then we define the dynamic thermal efficiency 

Et(Ts,l,Ps,l,Ps,l,{l) - E T (r R (l + Ag),p L ,xp L , 1) 



(13) 



(14) 



and the dynamic magnetic efficiency 

_ E M (r s _ L ,p sx , a SiL , grj - E„{T R {\ + Ag),p L , cr L , 1) 

S M,L ■- ~ , 

E 

where E is the total initial energy of both shells 
E := E K (r R (l + Ag),p L , 1) + E T (T R a + Ag),p L , X p L , 1) 

+E M (T R (l + Ag),p L ,cr L , 1) + E K (T R , 1, 1) 

+E T (T R ,l,xA) + E M (T R ,l,o- R ,l). 

Equations d 1 3 1 ) and d 1 41 > express the fraction of the initial en- 
ergy that the RS has converted into thermal and magnetic energy, 
respectively. 



(15) 
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3.2.2 Forward shock 

In complete analogy, we define the thermal and magnetic efficien- 
cies for the forward shock, 



S T,R 



£m,r 



E t (T s ,r,Ps.r, Psjt, fo) ~ E T (T R , l,x, 1) 
Em(Fs,r,Ps,r, o- S jt, £ R ) - E M (T R , 1, tr R , 1) 



(16) 



(17) 



where 



V S,R ~ V CD 



< 1 . 



Vsjt is the velocity of the FS in the lab frame. Here we set e tr = 
Bm.r = if the forward shock is absent. 

Combining equations d!3t . t!4h d!6t and (TTJ we define the 
dynamic thermal and magnetic efficiency of internal shocks 



e t :- e tl + b tr 
e m '■= e m,l + e m,r ■ 



(18) 
(19) 



We point out that these definitions of efficiency generalize the 
ones typically used when cold, unmagnetized shell collisions are 
considered. In that case, initially one only has bulk kinetic energy 
in the shells (i.e., Eq = E K (T L ,p L , l) + E K (T R , 1,1)). In case of colli- 
sions of arbitrarily magnetized shells with arbitrarily initial thermal 
content, E can be substantially larger than the initial kinetic energy 
in the shells. 



3.3 Efficiency of energy dissipation by a rarefaction 

In a rarefaction there is a net conversion of magnetic and/or ther- 
mal energy into kinetic energy, thus the net dynamic efficiency 
produced by a rarefaction, defined as in e.g., Eq. dl4t . should be 
negative. The consequence of this is that, when we obtain a shock- 
contact-rarefaction or rarefaction-contact-shock structure as a solu- 
tion to the Riemann problem, it may happen that the total (left plus 
right) thermal or magnetic efficiency (Eqs. d!8t-l|19t) was negative. 
However, this situation does not correctly model the fact that, in 
cases where a shock exists only in one of the shells, it is still able 
to radiate away part of the thermal or magnetic energy behind it, 
even though there is a rarefaction propagating through the other 
shell. We also point out that rarefactions happen in our case, where 
we model initially cold flows, as a result of a net conversion of 
magnetic into kinetic energy. This kinetic energy can be further re- 
cycled by the flow, and dissipated in the course of ulterior binary 
collisions. Therefore, directly summing the (positive) dynamic ef- 
ficiency of conversion of kinetic-to-thermal/magnetic energy in a 
shock to the (negative) dynamic efficiency of conversion in a rar- 
efaction is inadequate. The total dynamic efficiency in a case where 
only one shock forms will be determined only by the efficiency in 
the shocked shell. Thus, we set e TL = e m ,l = (e tr = e m .r = 0) 
if the reverse (forward) shock is absent. We note that this contrasts 
with previous works o n internal collisions of unmagnetized shells 
(e.g.. lKino et al 1l2004 . and may yield higher values of the net dy- 
namic efficiency. 



4 PARAMETRIC STUDY OF THE DYNAMIC 
DISSIPATION EFFICIENCY 

Next we study the dynamic dissipation efficiency in the process 
of collision of cold, magnetized shells.The shells are assumed to 




-6 -4-2 2 

logcr L 

Figure 1. Contours: total dynamic efficiency sj + em (eqs. U8t and il9> ) in 
the blazar regime (T R = 10, Ag = 1) for different combinations of (07,, <tr). 
Contours indicate the efficiency in percent and their levels are 1, 2, 3, 4, 5, 
6, 7, 8, 9, 10, 11, 12 and 13. In the region of the parameter space above 
the dashed line there is no forward shock, while the reverse shock is al- 
ways present for the considered parametrization. Filled contours: magnetic 
efficiency in percent. 



be cold because in the standard fireball model (e.g., |Piranll200^) 
, almost all the internal energy of the ejecta has been converted 
to kinetic energy before internal shocks start to show up. Thus, 
a regime where the parameter x ls large does not properly model 
an efficiently accelerated ejecta by non-magnetic processes. On the 
other hand, if the ejecta were accelerate d by magnet ic fields (like in 
Poynting-dominated flow models; e.g.. lUsovll 19921) . then the flow 
is cold all the way from the beginning to the internal shock phase, 
and then x should also be small in such a case. 

For all the models in this paper, in order to reduce the dimen- 
sions of the parameter space, we fix^ = 10~ 4 uniform everywhere, 
to model initially cold shells and, unless otherwise specified, set 
Ag = 1 as a reference value. We choose x sufficiently small so that 
it does not influence the solution of the Riemann problem. In the 
first two subsections we consider blazar and GRB regimes. Then, 
we study the flow structure for three representative Riemann prob- 
lems, and end the section with a discussion of the impact of varying 
Ag on the efficiency. 

4.1 Blazar regime 

In the blazar regime we set T R = 10 as a typical value of the Lorentz 
factor of the outflowing material. We continuously vary a L and a R 
and show contours of total efficiency (e t + e m ) in Fig.Q] 

The maximum efficiency is attained for moderately magne- 
tized slower shells (cr R as 0.2) and highly magnetized left shells 
(cr L x 1). The broad region to the right of the efficiency maximum 
is independent of cr L because in a collision with such a highly mag- 
netized fast shell almost all the energy is dissipated by the FS. In 
the region above the dashed line of Fig. [T] the FS is absent and, 
thus, since only the RS dissipates the initial energy, the efficiency 
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Figure 2. Same as Fig.[T] but in the GRB regime (T R = 100, Ag = 1). 



slightly drops. However, the transition between the regime where 
the two shocks operate or only the RS exists is smooth. The reason 
being that the efficiency below the separatrix of the two regimes 
and close to it is dominated by the contribution of RS. 

As expected, when either <r L or cr R approach low values, the 
dynamic efficiency ceases to depend on them. This can be seen in 
the center of the left side of Fig. [TJ where, for cr L <k 1, the dynamic 
efficiency only depends on cr R . The converse is true in the center of 
the lower side of the figure, where cr R <k 1 . 



4.2 GRB regime 

The results in the GRB regime (T R = 100, Ag = 1) are shown 
in Fig. [2] The general shape of the contours is similar to Fig. Q] 
which is expected since both in the blazar and in the GRB regime 
the flow is utrarelativistic and most of the quantities which de- 
pend on the relative velocity between the faster and the slower 
she ll depend only weakly on T R , Ag being the crucial parame- 
ter jDaigne & Mochkovitch][l998l) . For example, the dashed curve 
which delimits regions with and without a forward shock does not 
depend on T R but only on Ag. 

The maximum of the dynamic efficiency in the GRB regime is 
localized at roughly the same spot as in the blazar regime. However, 
compared to the former case the region of maximum efficiency (i.e., 
where b t + s M > 0.13) is smaller. 



4.3 Flow structure 

In this subsection we study the flow structure for three represen- 
tative models in the GRB regime. Their location in the param- 
eter space is marked by letters A, B and C in Fig. [2] Model A 
corresponds to a prototype of interaction between non-magnetized 
shells (cr L = cr R = 10~ 6 ). Model B is picked up to illustrate the 
flow structure at the maximum efficiency (cr L = 0.8, cr R = 0.2). 
Finally, model C corresponds to the case when the FS is absent 
(07, = 1, ctr = 10 2 ). We show the rest mass density profile of these 




Figure 3. Rest-mass density profile of models A, B and C (see legends) 
whose parameters are given in Sect. |4.3| Profiles have been shifted so that 
the CD for all models coincides. In models A and B the FS and the RS are 
clearly visible, while in the model C a rarefaction wave is visible as a small 
"step" to the right of the CD. 

models in Fig. [3] Model A (thick full line on Fig. [3} shows two 
strong shocks which dissipate kinetic into thermal energy. In con- 
trast, model B (dashed line) has much weaker shocks due to non- 
negligible magnetization in both shells. Finally, model C (thin full 
line) does not have forward shock due to very high magnetization 
in the slower shell. 

All three models have a substantial dynamic efficiency, but 
there is a qualitative difference among them. In model A internal 
shocks dissipate kinetic to thermal energy only (thermal efficiency). 
In model B shocks mainly compress the magnetic field (magnetic 
efficiency) and dissipate only a minor fraction of the initial kinetic 
and magnetic energy to thermal energy. Finally, in the model C only 
the reverse shock is active, compressing the faster shell magnetic 
field. 



4.4 Dependence on Ag and on As 

The choice of a relatively small value of Ag in the previous sec- 
tion is motivated by the results of numerical simulations of rel- 
ativistic outflows (e.glAlov et al.ll200d,]2004 IMizuta et al.l | 2004 
[ Zhang et al.ll2003l l2004lMorsonv et al.N2007l ; Lazzati et al]|2009t 
IMizuta & Aloy||2009l) where Ag < 2 between adjacent parts o f the 
flow that may catch up (but see, e.g., lKino & Takaharall2008l who 
find Ag ~ 1-19 appropriate to model Mrk 421). This adjacent flow 
regions can be assimilated to pairs of shells whose binary collision 
we are considering here. 

Ho wever, it has been confirmed by several independent works 
(KPS 97; lBeioborodovll2000l : iKobavashi & Sarill2000l : | Kino et al.l 
I2004L etc.) that, in order to achieve a high efficiency (more than 
a few percent) in internal collisions of unmagnetized shells, the ra- 
tio between the maximum (r max ) and the minimum (r min ) Lorentz 
factor of the distribution of initial shells should be r max /r m i n > 10. 

In view of these results, we have also made an extensive anal- 
ysis of the dependence of the dynamic efficiency on the variation of 
Ag. Since we are also interested in evaluating the influence of the 
magnetic fields on the results, we define a new variable 



As 



1 +cr L 
1 +cr R 



(20) 
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log Ag 

Figure 4. The gray scale indicates the value of the maximum total dynamic 
efficiency (in percent) as a function of the parameter pair (Ag. As). The val- 
ues of the rest of the parameters are fixed to T R = 100, and x = 10~ 4 . 
Contours: magnetization of the slowest shell: <t r = 0.1, 0.5, 1, 5 and 10. 

and plot (Fig. |4) the value of the maximum efficiency reached for 
every combination (Ag, As) and fixed values of the rest of the pa- 
rameters to T R = 100, and x = 10~ 4 . To be more precise, for 
fixed Ag and As we need to look for the maximum of the ef- 
ficiency of all models whose <t l and cr R satisfy Eq. J20t . It is 
evident from Fig. [4] that the maximum total dynamic efficiency 
grows (non-monotonically) with increasing Ag, in agreement with 
the above mentioned works (where unmagnetized collisions have 
been considered). Indeed, a large value Ag > 10 yields dynamic 
efficiency values ~ 40% if both shells are moderately magnetized 
(cr R ~ a L < 0.1). Nevertheless, the amount of increase of efficiency 
with Ag depends strongly on As. For |As| > 1, corresponding to 
cases where the slower shell is highly magnetized (cr R > 4), the 
maximum dynamic efficiency is almost independent of Ag; while 
for |A.s| < 1, the maximum dynamic efficiency displays a strong, 
non-monotonic dependence on Ag. 

It is remarkable that values 5 < As < 100 yield dynamic ef- 
ficiencies in excess of ~ 20%, regardless of the relative Lorentz 
factor between the two shells. In this region of the parameter space 
the maximal dynamic efficiency happens when both shells are mag- 
netized (cr R > 10, cr L > 50), and the total dynamic magnetic effi- 
ciency dominates the total dynamic efficiency. 



5 DISCUSSION 

We have focused in this paper on the estimation of the dynamic 
efficiency of conversion of kinetic-to-thermal/magnetic energy in 
collisions (internal shocks) of magnetized shells in relativistic out- 
flows. A fundamental difference between the internal collisions in 
magnetized and unmagnetized outflows is the fact that in the for- 
mer case not only shocks but also rarefactions can form. Thus, one 
would naturally expect a reduced dynamic efficiency in the magne- 
tized case. However, we have shown that such dynamic efficiency 



may reach values ~ 10% - 40%, in a wide range of the parame- 
ter space typical for relativistic outflows of astrophysical interest 
(blazars and GRBs). Thus, the dynamic efficiency of moderately 
magnetized shell interactions is larger than in the corresponding 
unmagnetized case. This is because when the shells are moderately 
magnetized, most of the initial shell kinetic energy is converted to 
magnetic energy, rather than to thermal energy. 

The difference in efficiency between flows with moderate 
Lorentz factors (T R = 10) and ultrarelativistic ones (T R = 100) 
is very small, because in the ultrarelativistic kinematic limit (i.e., 
T » 1), once the energy flux and the magnetizations of both shells 
are fixed, the key parameter governing the dynamic efficiency is Ag 
rather than T R . From numerical simulations one expects that any 
efficiently accelerated outflow will not display huge variations in 
the velocity between adjacent regions of flow. Therefore, values of 
Ag ^ 1 seem reasonable and Ag = 1 has been taken as a typi- 
cal value for both blazars and GRB jets. A fixed value of Ag = 1 
brings maximum efficiency when the magnetizations of the collid- 
ing shells are (<t l , <t r ) — (1, 0.2). Larger dynamic efficiency values 
~ 40% are reached by magnetized internal shocks if Ag > 10 and 
|As| <, 1, corresponding to cases where the magnetization of both 
shells is moderate (<r R a L < 0.1). 

Consistent with our previous work dMimica et alj|2007h . in 
the limit of low magnetization of both shells, the kinetic energy 
is mostly converted into thermal energy, where the increased mag- 
netic energy in the shocked plasma is only a minor contribution 
to the total dynamic efficiency, i.e., e t <k e m . Here we find that 
as the magnetization of the shells grows, the roles of e t and e m 
are exchanged, so that s T < e m (at the maximum dynamic effi- 
ciency e t ^ 0. \e m ). If the magnetization of both shells is large, 
the dynamic efficiency decreases again because producing shocks 
in highly magnetized media is very difficult. All these conclusions 
are independent on the EoS used to model the plasma, i.e., they 
are both qualitatively and quantitatively basically the same inde- 
pendent on whether a polytropic EoS with a fixed adiabatic index 
is taken (either y = 4/3 or y = 5/3) or a more general, analytic 
approximation to the exact relativistic EoS (the TM EoS; see Ap- 
pendix) is considered. 

The comparison of our results with previous analytic or semi- 
analy t ic works (e.g.,KPS97; | Beloborod ov 2000; Kobavashi & Saril 
1200 ll ; IKobavashi et"all |2002» ~ is not straightforward. Generally, 
these works compute the efficiency of the collision of shells with- 
out computing their (magneto-)hydrodynamic evolution and, on the 
other hand, these works include not only a single collision, but the 
multiple interactions of a number of dense shells. The bottom line 
in these previous works is that internal collisions of unmagnetized 
shells can be extremely efficient; the efficiency exceeding 40%, or 
even ~ 100% dBeloborodovfeOOOl) if the spread of the Lorentz fac- 
tor (i.e., the ratio between the Lorentz factor of the faster, r max , and 
of the slower r m j n shell in the samp le) is large (r max /r m i„ = 10 3 ; 
e.g., KPS97, IKobavashi et"aT]|2002l) . For a more moderate spread 
of the Lorentz factor r max /r m j„ = 10, the efficiency is ~ 20%. We 
note that these high efficiencies are reached because a large number 
of binary collisions is included in the model (not only a single one 
as in our case). Thus, the kinetic energy which is not dissipated in 
the first generation of collisions (between the initially set up shells), 
can be further converted into internal energy as subsequent gener- 
ations of collisions take place. In contrast, we find that moderate 
magnetizations of both shells (cr < 0.1) and Ag > 10 (which would 
roughly correspond to r max /r min = 9) are enough for a single bi- 
nary collision to reach a total dynamic efficiency of ~ 40%. 

We point out that the energy radiated in the collision of magne- 
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tized shells is only a fraction , f r ^0.1 (e.g ., Panaitesc tTet alj 19991 : 
iKumaJI 19991) to f r 1 (e.g.. lBeloborodovll2000r) of the energy dv- 
namically converted into thermal or magnetic energy. Thus, the ra- 
diative efficiency of the process, measured as the fraction of the to- 
tal initial energy converted into radiation, will be 1 //, times smaller 
than the computed dynamic efficiency. Even considering this factor, 
single binary collisions between moderately magnetized shells may 
yield efficiencies ~ 0.4 /j , which can obviously rise if many binary 
collisions take place in the flow reprocessing the remaining kinetic 
energy of the first generati on of interacting shells (i n the same sta- 
tistical way as discussed bv lKobavashi & Sar il200Ch . Therefore, on 
the light of our results, binary collisions in relativistic magnetized 
flows are efficient enough, from the dynamical point of view, to be a 
valid mechanism to dissipate the bulk kinetic energy of relativistic 
ejecta. Hence, the main restriction on the radiative efficiency comes 
from the radiation mechanism setting the limiting factor f r . 

We stress that we are not assuming any particular radi- 
ation mechanism in this study (thus, determining a value for 
f r ). Therefore, we compute the dynamic efficiency including not 
only the increased thermal energy in the flow, but also the ex- 
tra magnetic energy resulting from the magnetic field compres- 
sion. This is justified because, although standard shock accelera- 
tion mechanisms are ineffi cient in very magnetized shocks (e.g., 
ISironi&Spitkovskvll2009l) . other mechanisms might extract the 
energy from the whole volume of a very magnetized fluid (e.g., 
iThompsonll 19941 : iGiannios & Spruitll2007l) . 

The estimated dynamic efficiency in the binary collision of 
magnetized shells will be completed in a future work by account- 
ing for the numerical MHD evolution of such building blocks of the 
internal shock models. A step further would be to compute the ra - 
diative efficiency using the method devised in lMimica et al. 
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APPENDIX A: EQUATION OF STATE 

In this Appendix we discuss the effects of using a more re- 
alistic equation of state on our results. We choose the TM 
analytic approximation to the Synge equation of state (EoS) 
dde Berredo-Peixoto et af]|2005l : iMignone et al]|2005h . In the TM 
EoS the specific enthalpy can be written as (using the notation of 
SectionED 
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+ 1 



1/2 



and the specific internal energy 

K 2 -i 1/2 



2p 
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Figure Al. Same as Fig. [2] but for the TM equation of state. 

From J A 1 1 > it can be seen that the limit a = the effective adia- 
batic index of this E oS lies between 4/3 and 5/3. We modified the 
iRomero et all j2005h solver to include the TM EoS. 

On Fig. I Al I we show the dynamical efficiency in the GRB 
regime using the TM EoS. Comparison of Fig.|2]and Fig. lAll shows 
that, overall, the dynamical efficiency is higher when using the TM 
EoS, but the qualitative features remain the same in both cases. 
Also, as expected, in the highly magnetized regime the differences 
are minor, since in both cases (polytropic or TM EoS) the effective 
adiabatic index approaches 2 in such a regime. Figure IA2I shows 
the influence of the EoS on the existence of the FS. The only dif- 
ference between models with different EoS appears in the region 
where the faster shell is weakly magnetized. There a slightly higher 
(lower) magnetization of the slower shell is needed to suppress the 
FS when using TM than when using a polytropic EoS with y = 4/3 
(y = 5/3). The overlap of all three curves in Fig. lA2l in the limit 
of high magnetization of both shells, shows again the result that 
the choice of EoS plays no role in such a regime, as expected. We 
note that the separatrix between the regions of existence and non- 
existence of the FS corresponding to the case y = 5/3 lies closer 
to that of the TM EoS than that corresponding to the case y = 4/3. 
This is a natural consequence of the fact that the initial shells are 
both cold, and thus, the effective adiabatic index is closer to 5/3 
than to 4/3, at low magnetizations. 

We note that in the case of the external shocks, where rela- 
tivistic fluid encounters a cold external medium the choice of the 
realistic EoS can influence results muc h more dramatically than in 
the case of the internal shocks (see e.g.. lMeliani et alj2008t) . 
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